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Abstract 



We report on results of the effective theory method applied to neutron-deuteron 
scattering. We extend previous results in the J = 3/2 channel to non-zero energies 
and find very good agreement with experiment without any parameter fitting. 



Since the establishment of QCD as the theory of the strong interactions very little progress 
has been made in understanding nuclear forces from first principles. Many phenomenolog- 
ical models have been developed with great success, but they all suffer from shortcomings, 
among them ambiguities in using nucleon-nucleon scattering information in the calculation 
of other processes, difficulty in relating them to the underlying QCD, and especially, lack 
of a systematic expansion in a small parameter. The effective field theory approach has 
the promise of solving these difficulties [p]]. The role of the small parameter is played by 
the ratio of the typical momentum scale Q in the problem to the scale associated with the 
physics left out of the effective theory. In the case of nucleon interactions up to momenta 
of the order of 300 MeV one can build an effective theory containing only nucleons and 
pions (and delta isobars). The scale of the physics left out is ~ m p and the expansion pa- 
rameter is ~ Q/m p . This idea was elaborated in a large number of works in the last few 
years J2|. Subtle problems regarding the naturalness of the shallow nuclear bound states, 
renormalization, and power counting in the presence of pion exchange are nowadays sub- 
ject of intense discussion || |J. However, such problems can be bypassed in those nuclear 
processes where the typical momentum scale is small compared to the pion mass. In this 
case one is allowed to use an effective theory without explicit pions, contact forces being all 
that remains. That is what happens in deuteron physics, since the typical momentum scale 
in a deuteron is given by the inverse of the 3 Si scattering length, 1/at << m n . This situ- 
ation arises because the nuclear potential is fine-tuned so that there is a bound state close 
to threshold with energy ~ l/(Maf), much smaller than other scales in the problem like 
~ l/(MrQ t ) ~ m^/M (we take the effective range in the 3 Si channel r ot ~ l/m w for power 
counting purposes). Attempts at model- independent approaches in this energy range have a 
long history ||. When this approach is applied to nucleon-nucleon scattering up to momenta 
~ l/a t the effective range expansion is reproduced. The first non-trivial application is thus 
in the three-nucleon sector. In this rapid communication we report results of this approach 
in the case of neutron-deuteron scattering in the J = 3/2 channel below deuteron break-up. 
We perform an expansion on powers of rot/cit and pro t , where p is the typical momentum 
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of the process, keeping terms up to order (r 0t /a t ) 2 , (pr 0t ) 2 (we take r ot ~ l/m n ). Results in 
extraordinary agreement for the quartet scattering length were previously reported in Ref. 
||. Here we extend this calculation to finite energy. 

In the J = 3/2,1 = 1/2 channel the spins of all three nucleons are aligned and all two- 
body s-wave interactions are in the spin triplet, isospin singlet channel. (For this reason we 
will drop from now on the subscript in 04 and rot). The effective Lagrangian restricted to 
this channel is given by 



C = N \id Q + ^ + ...)N + C (N^T 2 aa 2 N) 2 (1) 



+ C 2 \(N ] r 2 aa 2 VN){N ] r 2 aa 2 VN) - 3{N ] r 2 aa 2 N)(N ] T 2 aa 2 V 2 N) + h.c. 



+ ... 



where M is the nucleon mass, C n are constants related to the two-body force terms containing 
n derivatives, and the dots stand for higher-order terms including relativistic corrections, 
higher- derivative terms, three-body forces, etc. The constants C n are determined by nucleon- 
nucleon scattering data. It turns out that, using dimensional regularization and minimal 
subtraction, C ~ a/M, C 2 ~ r (r a) / M , C 4 ~ r (r a) 2 / M + . . . and so on (ellipses stand 
for terms suppressed by powers of r /a ). The leading pieces in each one of these terms form 
a geometric series that can be conveniently summed to all orders by the introduction of a 
field of baryon-number two |7J 



£ = N ^' d o + ^ + ---)N + ^-Hdo-J^ + A + ...)d 

■Naa 2 T 2 N + h.c.) + ... (2) 

More generally, if the dibaryon field d is integrated out, the Lagrangian (§) is recovered as 
long as A and g are appropriate functions of Co and C 2 . This resummation is by no means 
necessary, since for momenta of the order p ~ 1/a the resummed terms are subleading, but 
it is a convenient way of computing higher-order corrections. 
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Figure 1: Dressed dibaryon propagator. 

The numerical values of g and A can be determined if we consider the dressed dibaryon 
propagator (Fig. p. The linearly divergent loop integral is set to zero in dimensional 
regularization and the result is 

iS (p) = ~ 2 M a ) • ( 3 ) 

This propagator is, up to a constant, the scattering matrix of two nucleons in the 3 Si channel, 

T(yk ^ = ~M 2nA 2n P i}? ^ 

where k 2 /M is the energy in the center-of-mass frame. This result is just the familiar effective 
range expansion, from what we can infer the proper values for the constants g and A. Using 
a = 5.42 fm and r = 1.75 fm 0, we find 

g 2 =J fa =1.6.10- s MeV-\ (5) 
^ =UW =8-7MeV. (6) 

From Eqs. (|3]), (J5|), and @ we see why it is necessary to resum the bubble graphs in Fig. [I] 
to all orders for p ~ 1/a: the term in the square root coming from the unitarity cut is of the 
same order as A. On the other hand, as mentioned before, the kinetic term of the dibaryon 
is smaller than the other terms in @ and is resummed for convenience only. Notice that 
the propagator (|3]) has two poles, one at p° = p 2 /4M — B (the deuteron pole), another at 
p° = p 2 /AM — Bdeep (unphysical deep pole), and a cut along the positive real axis starting 
at p° = p 2 /AM. 

Let us now turn to neutron-deuteron scattering. The simplest diagram contributing to this 
process is the first diagram in Fig. 0. For momenta of the order of p ~ 1/a it contributes 
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Figure 2: Graphs contributing up to order (r /a) 2 . 

~ Mg 2 /p 2 ~ a 2 /Mr . The one-loop graph mixes different orders of the expansion, since 
it involves the dibaryon propagator g 2 /{A + p 2 /M) ~ (o/M)(l + O(ro/a) + ...); it gives a 
contribution ~ g 4 M 2 /pA ~ (a 2 /Mro)(l + O(ro/a) + ...). It is easy to see that the remaining 
graphs in Fig. ^ give contributions of the same order, which means that an infinite number 
of diagrams contribute to the leading orders. 

Other contributions are suppressed by at least three powers of r /a or pr |3J. For instance, 
the effect of the subleading (not resummed) piece of C4 is to generate the shape parameter 
(~ k 4 ) in the effective range expansion of the nucleon-nucleon interaction. Its typical size is 
~ k 4 r$ compared to the leading piece ~ 1/a and is thus also suppressed by (r /a) 3 . Likewise, 
p-wave interactions, unaffected by the existence of a shallow s-wave bound state, arise from 
a term in the Lagrangian with two derivatives and a coefficient of the order ~ l/MmJ. 
We conclude then that a diagram made out of the substitution of one of the dibaryon 
propagators in a diagram in Fig. |] by a p-wave interaction vertex would be suppressed by 
(ro/a) 3 in comparison to the leading order. Three-body force terms have to contain at least 
two derivatives since in the J = 3/2 channel all the spins are up and Fermi statistics forbids 
the placement of all three nucleon in a s-wave. The natural size of the coefficient of the six 
nucleon, two derivative term that produces such a three body force is ~ 1/MmJ. This term 
is generated, upon integration of the dibaryon field, by a term containing two dibaryon fields, 
two nucleon fields and two derivatives with a coefficient of the order of ~ r\jMa 4 . Thus 
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contributions coming from the three-body force are suppressed in relation to the leading 
order graphs by (r /a) 6 . 

A calculation accurate up to corrections of order (ro/a) 3 is possible by summing the dia- 
grams of Fig. |2|. Fortunately, the interaction mediated by the s-channel dibaryon generates 
a very simple, local and separable potential between nucleons. It is well known that the 
three-body problem with separable two-body interactions reduces to an equivalent two-body 
problem. In our case the equation to be solved can be read off Fig. 0, and an integration 
over the energy inside the loop gives || 
3(p 2 - k 2 ) 



8MV 



1 

47T 



-1 



-(p 2 - k 2 ) + MB - VMB) 



dH 



tjp, k) 
-k 2 - 



(7) 



IE 



t(l,k) 



(p- k/2) 2 + MB J (^fl 2 -l-p + p 2 -\k 2 + MBl 2 -p 2 -ie 
where B is the deuteron binding energy. Since we are interested only in s-wave scattering, 

we should project this equation into its L = component. The result is 



-T] + 

2 

TlX JO 



1 



Ux 2 -y 2 ) + l 



a(x, y) = In 

xy 



dz zln 



' z 2 + x 2 + 1 — \y 2 + xz" 



z 2 + x 2 + 1 



\x + y/2) 2 + T 
(x-y/2) 2 + y 

a{z,y) 



3„,2 



y* — xz l z' 



y 



ie 



1 1 4MB, and 



where we use the dimensionless quantities x = p/VMB, y = k/vMB, z 
a(x,y) = ^¥ B tL=o(p, k), and rj = y/ MBr /2. For finite values of k this equation is complex 
even below threshold (3k 2 /A = B) due to the ie prescription. It is convenient for the 



numerical treatment to use the real i^-matrix defined by 

a(x,y) 



K(x,y) 



1 + iya(y,y)' 



which satisfies the equation 

K(x,y) = - 

with 

h(x,y,z) 



2 r°° 

h( x ,y,y) — / 

7T JO 



V 



1 



xzf(x,y) 



-In 



dz z 2 h(x, y, z) 



' z 2 + x 2 + 1 - \y 2 + xz" 



;K(z,y), 



(9) 



(10) 



z 2 + x 2 + 1 



\y 2 



xz 



The phase shifts can be obtained directly form the on-shell i^-matrix : 



kcot5 = - ff^-r — . (12) 
K( k k ) 

Defining f(x,y) by the equation 

f( x h{x,y,y) 2 [°° 2 f h{x,y,y) \ f(z,y) 

f{x, y) = T7 7 - - / dzz h(x, y, z) - — -%, y, z) — -, 13 

h(y,y,y) it Jo \ h{y,y,y) )z 2 -y 2 

the on-shell fT-matrix can be obtained by 

K(y,y) =-h(y,y,y) (14) 
+ - f dz(z 2 h(y,y,z)f(z,y) -y 2 h(y,y,y)f{y,y)) ) . 

y TV JO ^ ' Z A — y z J 

Rewriting Eq. (P) this way greatly simplifies its numerical solution, for now the integrand 
is regular and the principal value can be dropped from Eqs. ( |i~3"l) and fli~4l) . 

We have solved Eqs. ([13]) and (|14|) numerically and the result for the phase shifts for 
energies up to the break-up point is shown in Fig. |3|. The data points at finite energy were 
taken from the phase shift analysis in and the much more precise (nearly) zero-energy 
point from ]ID) . Also plotted is the result of the leading order calculation obtained by setting 
r] = 0, in which case our equations reduce to the case studied in Ref. |T1| . 

We expect errors in our calculation to be of the order (r /a) 3 , (kr ) 3 compared to the 
leading order. These errors are smaller than the experimental errors in the finite energy 
case and of the same order as the experimental uncertainty in the case of the more precise 
measurement near k — 0, where we find 4 a t h = 6.33±0.10 fm H compared to 4 a = 6.35±0.02 
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Our results seem to deviate from a simple effective range type expansion only around the 
pole at ~ 0.05 fm -2 . (A pole in k cot<5 corresponds to a zero in the scattering matrix, 
which does not carry any special meaning.) This pole does not appear in potential model 
calculations (e.g., ||12||), and presumably will be smoothed out by higher-order terms that 
we have not yet included. It is interesting that the only "experimental" point in this region 
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Figure 3: k cotS in the J = 3/2 channel to order (r /a)° (dashed line) and (r /a) 2 (solid 
line) . Circles are from the phase shift analysis in Ref . and the triangle is from Ref . [[UJ . 

seems to indicate some structure there, but more experimental input would be necessary to 
confirm the behavior we predict. 

The calculation of higher-order corrections involves the knowledge of further counterterms 
like the ones giving rise to p-wave interactions, etc. These parameters can be determined 
either by fitting other experimental data or by matching with another effective theory - 
involving explicit pions — valid up to higher energies. If more precise experimental data - 
particularly at zero-energy — appear, we would be facing a unique situation where precision 



calculations in strong-interaction physics can be carried out QT3 1 and tested 
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